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As shown by Teukolsky, the master equation governing the propagation of weak radiation in a 
black hole spacetime can be separated into four ordinary differential equations, one for each space- 
time coordinate. ("Weak" means the radiation's amplitude is small enough that its own gravitation 
may be neglected.) Unfortunately, it is difficult to accurately compute solutions to the separated 
radial equation (the Teukolsky equation), particularly in a numerical implementation. The funda- 
mental reason for this is that the Teukolsky equation's potentials are long ranged. For non-spinning 
black holes, one can get around this difficulty by applying transformations which relate the Teukol- 
sky solution to solutions of the Regge- Wheeler equation, which has a short-ranged potential. A 
^— ^ ' particularly attractive generalization of this approach to spinning black holes for gravitational radi- 

■ ation (spin weight s = —2) was given by Sasaki and Nakamura. In this paper, I generalize Sasaki and 

' Nakamura's results to encompass radiation fields of arbitrary integer spin weight, and give results 

. directly applicable to scalar (s = 0) and electromagnetic (s = —1) radiation. These results may 

; I ■ be of interest for studies of astrophysical radiation processes near black holes, and of programs to 

compute radiation reaction forces in curved spacetime. 



o 
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I. INTRODUCTION 



(N 
> 

iH- ; In 1973, Teukolsky M derived a single partial differential equation describing the evolution of perturbations to 
rotating (Kerr) black holes. This master equation gives the linearized evolution of fields that arise from a perturbing 
source of stress energy — the charge and current densities associated with the perturbation — to the (vacuum) black 
, hole background. The solutions of the homogeneous version of this equation describe the propagation of radiation 
in black hole spacetimes. Thus, a common use of this formalism is to study the radiation emitted by matter in the 
environment of a black hole. In some cases, one can use such an analysis to study back reaction, determining how the 
O ■ perturbing source evolves as radiation carries away energy and angular momentum. A beautiful feature of the master 
equation is that it describes radiation fields of arbitrary spin weight s. It has been used extensively to study scalar 
5-H (s = 0), electromagnetic (s = ±1), and gravitational (s — ±2) radiation in Kerr spacetimes. 

The master equation is often|^ solved by introducing a multipolar decomposition of the radiation field. The solution 
^ I separates into functions of the Boyer-Lindquist coordinates: 

X ; = E Rimu.{r)sSZ{0y"'^e-'^' ; (i-i) 

each function is governed by an ordinary differential equation. (The precise meaning of is described in Sec. Q.) 
The t and (j) dependences are trivial. The 9 dependence is more involved, but can be evaluated in a straightforward 
matter. The functions sSf^ are spin-weighted spheroidal harmonics, which are generalizations of spin-weighted 
spherical harmonics to a spheroidal geometry. The spin-weighted spherical harmonics in turn are generalizations 
of spherical harmonics that encode the rotation properties of spin s fields; see Refs. for further discussion. A 
detailed algorithm for computing sSZ{0) is given in Ref. g. For the purposes of this paper, the 9 dependence is 
considered known. 

The radial dependence, Rimu:{r) = R{r), on the other hand, can be rather difficult to calculate in practice, particu- 
larly in a numerical computation. The fundamental reason for this difficulty is the nature of the equation that governs 
R(r): this equation (the Teukolsky equation) has a long-ranged potential. In source-free form, it can be written 

^ + FT{r)^+[^'-Ur{r)]R = 0, (1.2) 



^Separation is not always used. There is also a body of work that uses the master equation to evolve initial data. This 
approach has been extensively used to study the endpoint of binary black hole collisions; see Ref. jij and references therein. 
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where lo is the frequency of the radiation mode and 



2Mr+ r-r+ 2A-fr_ r-r^ 
T = r H in—— in——— (1.3) 

is the Kerr "tortoise coordinate". The potentials UTir) and Fxir) are rather complicated; they encode the most 
interesting features of wave propagation in black hole spacetimes, such as scatter from spacetime curvature^ and 
superradiant scattering (radiation whose scattered amplitude exceeds the ingoing amplitude due to extraction of 
energy from the black hole's spin). For large r, 

, , AisLU X + 2amuj + 81 M suj o 
Urir)^ — + -2 + 0{l/r^). (1.4) 

[The quantity A is related to the eigenvalues of the 9 dependence; see Ref. ||] for details.] For l arge r, Fx^r) and 
Uxir) fall off only as l/r — they are long-ranged, like the Coulomb potential. The solution of Eq. ( |l.2| ) for large r is 



R-C,'-^+C2^. (1.5) 



The complex constants Ci and C2 are determined by boundary conditions. This asymptotic solution illustrates the 
difficulty in solving the Teukolsky equation: the coefficient of e*'^'' differs from the coefficient of e"*'^'' by r^^*. For 
negative s, this becomes extremely large — large enough that the ingoing e~^'^^ piece will eventually be entirely 
lost in any numerical computation due to round-off error. Hence, for negative s, it is nearly impossible to set proper 
boundary conditions on the solution's phase at large r. Similarly, for positive s it is difficult to set boundary conditions 
near the event horizon^. For r very close to = M + \J — a?- (the location of the event horizon in Boyer-Lindquist 
coordinates), the solution is 

R = CsA-^e-^P"' -f C4e'f''* . (1.6) 

I have introduced p = w — muj^^ where — a/2Mr^ is the angular velocity at which observers at the horizon are 
seen to rotate. Because the Boyer-Lindquist coordinates t and (j) become twisted and entangled near the horizon, p 
describes the frequency of wave modes in that region. The factor A = — 2Mr + a? goes to zero at the event horizon. 
Hence, for positive s, the ingoing solution swamps the outgoing solution as one approaches the horizon. Whether one 
chooses positive or negative s, there exists a domain in which one cannot accurately compute numerical solutions by 
directly integrating the homogeneous Teukolsky equation. 
Various approaches have been discussed to circumvent this difficulty. One of the firs t wa s introduced by Teukolsky 



and Press ||7|,§|- Their approach used the fact that, for a given |s|, the solutions to Eq. (1^) for s = +\s\ and s = — |s| 
are physically equivalent: there exist rules to take the positive s solution to the negative s solution, and vice versa. 
Thus, one can initially choose s so that the solution is well-behaved in the initial r domain, and then "switch horses" 
and integrate with the other sign of s as the integration approaches the other asymptotic domain. 

A somewhat more elegant way to compute R was developed by Chandrasekhar B. As already noted, the poor 



behavior of the solutions (1.5) and ( |1.6| ) is due to the long-ranged nature of the Teukolsky equation's potentials. 



Rather than try to work with an equation that is simply not well-behaved, one should find transformations which 
relate the Teukolsky solution R to the solution X of some equation whose potentials are short-ranged. For example, 
when the black hole spin is zero, black hole perturbations can be described using the generalized Regge- Wheeler 
equation |10[: 



dr 



,2 < Lu'-VRwir,s) X = 0, (1.7) 



■^As shown by Leonard and Poisson the phenomenon of tails (delayed propagation due to scatter from spacetime curvature) 
is to leading order independent of s, and is encoded in the logarithmic behavior of r* , not tire potentials. 

^One can expand the potentials near the horizon and see that they die away slowly as r* —00 (which corresponds to 
r r+). The actual form is somewhat messy, and is not given here explicitly. 
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where 

..... 1(1 + 1) 2(^2 - 1)M' 
VRw[.r,s) = f 



(1.8) 



(Here, / = 1 — 2M/r.) This potential dies aw ay fa ster than l/r, and so is sh ort ranged. Chandrasekhar showed 



(for specific choices of s) that solutions to Eq. (1.7) and solutions to Eq. (1.2) (with a = 0) are related by simple 
rules. (Below I generalize these rules to any value of s.) This is extremely useful for numerical work: one can 
integrate Eq. ( |l.7| ) to accurately compute X and then transform to R. Note that in Ref. [|| Press and Teukolsky 
introduced a transformation that, in essence, transformed to a function governed by an equation with a better behaved 
potential. They did not, however, discuss the nature of the transformation in terms of the rangedness of the potentials; 
Chandrasekhar appears to have been the first to systematically approach this problem with the viewpoint that the 
long-ranged potential was the key issue. Note also that Chandrasekhar's notation H is rather different from that 



used here; I use a notation similar to that used in [y_l|. The transformation given in Refs. y] and |11| is for s = —2; 
a rule for s = — 1 is given in |^ . 

For spinning black holes, perhaps the most elegant generalization of Chandrasekhar's approach was given by Sasaki 
and Nakamura [Q. They derive a transformation rule which relates R for any physical spin a to the solution X of 
an equation whose potentials are short ranged. The transformation and short-ranged potentials are designed such 
that if a = 0, the potentials reduce to the Regge- Wheeler potentials. This approach is very natural in the sense that 
its solutions are monotonic with respect to spin, ranging from the Schwarzschild value a = to the extreme Kerr 
limit a = M. Chandrasekhar and Detweiler also investigated several transformations relating i? to a short-ranged 



solution X [13-16|. In some (but not all) cases, these transformations are governed by equations which reduce to the 
Regge- Wheeler equation in the a = limit; however, the equations themselves are often not as "nice" to work with. 
For example, the equation for X is sometimes given in terms of a f requ ency dependent variable f*{Lu) (see Ref. [|5]) 
which is different from the "usual" tortoise coordinate r* [cf. Eq. ( |l.3| )]. This coordinate can be doubly valued and 
mask features such as superradiant scattering. Also, the potentials of each of their equations are pathological for some 
set of frequencies |l7|. By using a set of multiple perturbation equations and transformation laws, one can always 
find a non-pathological tool for any given frequency. But, there is no single rule that works for all frequencies. (These 
difficulties do not mean that Chandrasekhar and Detweiler's approaches are not useful. Campanelli and Lousto |lj] 
used rules very similar to Chandrasekhar and Detweiler's in order to show that solutions to the Teukolsky equation 
are well-behaved even for sources that extend to infinity.) 

Sasaki and Nakamura's work is restricted to the choice s = —2. This is an appropriate choice for studies of gravita- 
tional perturbations, and so the Sasaki-Nakamura equation has been extensively used in studies of gravitational-wave 
generation and gravitational radiation reaction ||^,|l^,^-22 . Other values of s are interesting as well. For example. 



s = ±1 corresponds to electromagnetic radiation. The propagation and scatter of electromagnetic waves in black hole 
spacetimes is of great astrophysical interest. Also, much work is currently being directed toward understanding how 
one calculates self forces and radiation reaction forces in curved spacetimes |23-3^]. Implementations of the general 



formalism (Refs. ||2^,|2J]) to date have been restricted to scalar (s = 0) or electromagnetic fields. They have also been 
restricted to spherically symmetric spacetimes. Tools for effective calculation of radiation fields in Kerr spacetimes 
will make it possible to extend these calculations to more realistic spinning black holes. 

In this paper, I generalize the Sasaki-Nakamura equation to arbitrary integer spin weight s. The generalized 
Sasaki-Nakamura (GSN) potentials are given for any s, but in terms of two unknown functions, a{r) and /3(r). These 
functions are fixed by requiring that the transformation which relates the Teukolsky solution to the GSN solution 
reduces, in the Schwarzschild limit, to the transformation between the Teukolsky solution and the Regge- Wheeler 
solution. They also must be chosen so that the potentials they generate are of short range. 

There is a great deal of freedom in how one chooses a and /3: for each value of s, there are an infinite number of 
functions which lead to a transformation with the correct limiting value and that produces short-ranged potentials. 
It is thus most practical to develop a and /3 on a case by case basis, rather than trying to develop generic formulas. 
Given the interest in scalar and electromagnetic radiation, I provide examples of a and /3 for s — and s = —1. 

Throughout this paper, a prime denotes d/ dr, where r is either the Boyer-Lindquist or the Schwarzschild coordinate 
(which coordinate should be clear from context). An overbar denotes complex conjugation. The function A = 

— 2Mr -f- a^, and / = 1 — 2M/r. Thus, for Schwarzschild holes, A — r^f. The function R will always denote the 
solution to the homogeneous Teukolsky equation; X will always refer to the solution of the equation with short-ranged 
potential. Section ^ reviews important aspects of the Teukolsky equation and its solutions, particularly how one can 
compute the solution given an appropriate source term and solutions to the homogeneous equation, and how those 
solutions are related to physical radiation fields. In Sec. p^ , I review the transformation rules for Schwarzschild 
black holes. These rules, which are particularly simple, will be used as guidelines for constructing transformations 



appropriate to radiation in Kerr spacetimes. Finally, in Sec. IV I construct the short-ranged equation for Kerr black 
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holes and provide a recipe for specifying the transformation rule for a radiation field of arbitrary spin weight. I apply 
this recipe in Sees. and VI to scalar and electromagnetic radiation fields, respectively. The resultant transformation 
rules and equations shoul d for m a useful basis for further studies of radiation in Kerr spacetimes. Some concluding 
discussion is given in Sec. 



VII 



II. SOME PROPERTIES OF THE TEUKOLSKY EQUATION AND ITS SOLUTIONS 



As background for the calculations in this paper, I review in this section the most important properties of the 
Teukolsky equation and its solutions. 

As discussed in the Introduction, Teukolsky 0] showed that one can separate the wave equation for a field of 
s pin weight s radiation propagating on a Kerr black-hole background using the multipolar decomposition given in Eq. 
(1.1). The 4> and t dependence is trivial, and the 9 dependence is straightforwardly dealt with. The r dependence, on 
the other hand, can cause problems. 



The radial function R{r) is a solution to the Teukolsky equation, Eq. (1.2). Here I write the Teukolsky equation 
with its source term and in terms of derivatives with respect to r rather than r*: 



(A'+^R')' - VTir)R = -r{r) . 

This is the way the Teukolsky equation usually appears in the literature. The potential Vt(?') is 

K{rf ~ 2is{r - M)K{r) 
VT{f) = AiSLor ; 



(2.1) 



(2.2) 



the quantity A = Sim ~ 2amuj + a}uj^ — s(s + 1), where Eim is the eigenvalue of the spheroidal harmonic [see in 
the Schwarzschild limit, = + !)]■ The function K(r) = (r^ + a^)a; — ma. 

The source term T(r) depends upon the spin weight of the radiation. It is constructed by projecting the radiation 
source onto legs of the Newman-Penrose null tetrad, 1, n, m, and m. A useful representation of the tetrad in 
Boyer-Lindquist coordinates is |1m 



1, — — , 0, —a sin 9 



A 



1,0,- 



a A sin^ 6'"' 



V2(r + ia cos Q) 



[iasin6',0, -E,-i(r^ 



2 , 

a sm 



in Q\ 



(2.3) 



The tetrad legs 1 and n represent ingoing and outgoing null vectors, respectively. Quantities constructed by project- 
ing onto 1 correspond to ingoing radiation and their sources; they map to positive s. Likewise, projection onto n 
corresponds to outgoing radiation and their sources^, and map to negative s. See Ref. for details. 

Because Eq. (2.1) is in self-adjoint form, one can construct its solution by the method of Green's functions 
This means that one needs to know only the solutions to the homogeneous equation. 



A- 



(A^+ii?')'-VT(r)i? = 0, 



(2.4) 



in addition to the source. One does this by adapting the generic solution, given in Eqs. (1.5) and (|1.6|), to the 
appropriate boundary conditions: no radiation may come in from infinity and none may come out from the event 
horizon. In other words, there exist two solutions, R^ (r) and R°°{r), whose asymptotic forms are 



R"{r) = 5*^°"= A 



hole A —s 



= B° 



r 



2s+l 



B" 



OO 



(2.5) 



(2.6) 



^Because one can transform between positive and negative s solutions, one can actually develop both ingoing and and outgoing 
radiation with a single source term. 
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The solution to the inhomogeneous equation ( |2.l[ ) which one constructs from Eqs. (p.5|), (2.6), and the source T(r) is 
conveniently written 

i?(r) = Z"{r)R°°{r) + Z°°{r)R"{r) , (2.7) 

where 



Z"{r)^—-^——f dr' A{r'YR"{r')T{r') 



Z°°{r) = ^ ^. ^ I dr' A{r'yR'^{r')T{r') . (2.8) 



Using Eq. (2.7) one can construct g^/. This quantity is related to a radiation field of spin weight s; the details of 
that relation depend upon the value of s. Typically, s'i' is constructed by projecting a tensor describing the radiation 
onto legs of the Newman-Penrose null tetrad. For example, o'l' = <i>, a massless scalar field. No projections are needed 
in this case. For s = ±1, we have 

-i^- = (r - ia cos0f (1)2 = (r - ia cos 9 f F^.^n^'m" , (2.9) 

where F^j/ is the electromagnetic field tensor. [There is a third projection, 0i = (1/2)F^^(P?t,'^ + fh^^my). It does not 
describe the radiative degrees of freedom of the electromagnetic field, and so is of less interest here.] For s — ±2, the 
radiative quantities are 

-2* ={r - iacosOf^i = -{r - iacosefCap-y5n°'m'^n"'m^ . (2.10) 

The tensor CapjS is the Weyl component of the spacetime's curvature. The quantities ■0^, with i an integer from to 
4, are the Newman-Penrose projections of the Weyl curvature (see Ref. [^). For unperturbed black hole spacetimes, 
all components except %p2 = —Cafj-ysl'^rnPTri^n^ = M/(r — iacos9)^ can be set to zero with an appropriate choice of 
gauge. This is the non-radiative "background" component of the curvature; the perturbations "00 ^^nd "04 represent 
radiation on the background. 

The solution for the (linear) evolution of radiation of spin weight s in a Kerr black hole spacetime is thus completely 
described by construction of the source T(r) appropriate to that spin weight and construction of the homogeneous 
solutions R^ (r) and R°°{r). As discussed in the Introduction — and, as should be clear from the asymptotic solutions 



( p.5[ ) and (2.6) — it is very difficult to build these solutions in a numerical integration. The remainder of this paper is 
devoted to methods for constructing R^ (r) and R°° (r) by finding transformations that relate the Teukolsky solution 
R{r) to solutions of equations with short-ranged potentials. 

III. RESULTS FOR SCHWARZSCHILD HOLES 

The Teukolsky equation for Schwarzschild black holes is 

{A'+^R')' ^VTs{r)R^O , (3.1) 

where 

Vrsir) = A - Aisruj + [2i.s{r - M)lu - {rLof]/f , (3.2) 
and A = A(a = 0) = l{l + 1) - s{s + 1). 



We would like to find rules that allow us to obtain R given a solution X of the Regge- Wheeler equation, (1.7). To 
do so, first define the quantity 

X - ^ . (3.3) 
If X satisfies the Regge- Wheeler equation, it is straightforward to show that x satisfies 
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(3.4) 



where 



By direct substitution, one can show that % can be transformed to i?, and vice versa, via 



(3.5) 



s < 

s = 
s > 



i? 



= ( r\/A I 



R 



^1 



X 
R 



R 



- I V 

r 



where 



X>± ^d/dr±iuj/J. 



(3.6a) 
(3.6b) 

(3.6c) 
(3.7) 



For s = — 1 and s = —2 , Eg. (3J)a) reduces to the Chandrasekhar transformation (see Ref. for s = -2, Ref. § 
for s = —1). Equations ( 3.6a ) - ( 3.6c ) serve as guidehnes that will be used to fix the form of the transformation rules 
for Kerr black holes. Note that the transformations from i? to x can be written 



(3.8) 



by repeatedly using Eq. (^^) to eliminate derivatives of second order and higher. The resulting functions a and [3 
may become rather complicated, particularly for large values of |s|, but the general operation is straightforward. (The 
factor A'^^^ is inserted for later convenience.) 



IV. PERTURBATION EQUATION FOR KERR HOLES 



Guided by Eq. (B.8), let us assume that functions a and (3 can be f ound that transform the Kerr solution R to 



solutions X of some other equation. By generalizing the relation (3.2) to a form appropriate for the Kerr metric 
and rewriting all derivatives in terms of r* we will come to an equation with short-ranged potentials governing the 
behavior of a function X{r). This is the generalized Sasaki-Nakamura (GSN) equation. It will depend explicitly on 
the (currently unspecified) functions a and B. T hese fun ctions will be specified by requiring that the transformation 
rule satisfy a form which reduces to Eqs. ( 3.6a ) ~ ( 3.6c ) when a — Q. This guarantees that solutions to the GSN 
equation are equivalent to solutions of th e R cgge- Wheeler equation in the Schwarzschild limit. 

To begin, differentiate x s-^id use Eq. (2.4) to eliminate the second derivative of R. The resulting equations for x 
and x' can be gathered neatly into matrix form: 



X\_( a /3A^+i \(R 

X' \a'+ fiVrl^' a + /S'A'^+i \R' 



A nice feature of Eq. (fl.lf) is that the inverse solution is rather obvious: 



If a + /3'A^+i 
ri (^-(a' + zJ^rA" 



-/3A^+i 
a 



(4.1) 



(4.2) 



where 



r; = a (a + /3'A"+i) - {a' + fiVrl^') 



(4.3) 
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is the determinant of the matrix in Eq. (4.1). 

Differentiating again and massaging the resultant expression gives us a second-order differential equation for x- 



A-^ (A^+^X')' - A.Fi(r)x' - .C/i(r)x = 
The potentials sFi{r) and sUi{r) are given by 

sFi{r) = -q'/i] , 



(4.4) 



,Uiir)=VT + 



/3A2 



(2a + /3'A^+i)'-^(a + /3'A«+i) 



Next, generalize Eq. (3^) to the Kerr form 



X 



X 



v/(r2 +a2)A* 



(4.5) 



(4.6) 



Using this to replace x for ^ in Eq. (4.4) and then replacing derivatives in r with derivatives in r* with the rule 

d _ jr^+a^) d 

dr A dr* ^ ' ^ 

yields the GSN equation: 

(4.8) 



The potentials are 



sFir) 
sU{r) 



dr* 



AsUijr) 
(r2 + a2)2 



d'^X , dX 
—^-sFir)—~Mir)X^O 



dr* 



The function sG{r) is 



Mr) 



2 AdsG/dr AMr)sFi{r) 

■^^y'^' + „2 , „2 „2 , „2 



rA 



s{r- M) 
(r2 + a2)2 ^ 7-2 + a2 



(4.9) 



(4.10) 



the functions sFi{r) and s?^i(?') are from Eq. (4.5). When s = —2, all functions reduce to those given by Sasaki and 
Nakamura (see Ref. [0). 

All of the quantities which have been derived to this point depend upon the as-yet-undetermined functions a and 
(3. We fi x the se funct ions by requiring that they affect a transformation between R and x which, as a — s- 0, reduces 
to Eqs. ( 3.6a ) - ( 3.6c ). A useful generalization of these transformations is 



s < : 
s = : 
s > : X 

where the operator 



X = (V(?-'+a')A) 9o{r)J- 

X = 9oir)R , 

-r 9o{r)J+ 



9i{r)J- 



92{r) ...J- 



s\{r)R 



9iir)J-^ 



92{r) . . . J+ 



9s{r) 



A 



yr^ + a^ 



R 



J± = d/dr±iK{r)/A 
generalizes 2?± to Kerr. The s — +2 transformation rule, for example, is 



X = 5o(0 J-f 



9i{r)J^ 



92{r)- 



:R 



(4.11a) 
(4.11b) 
(4.11c) 

(4.12) 
(4.13) 
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an example for s = — 1 is given in Sec. V] 



To now specify a and /3, one mus t pick fu nction s gi{r) and then repeatedly use Eq. (2^) to eliminate derivatives 
of second order and higher in Eqs. (4.11a) - ( [lllc| ). The resultant expressions for a and (3 will be, in general, quite 
complicated; examples are discussed in Sees, [v] and VI. The functions gi{r) must be chosen so that they become 
constant in the Schwarzschild limit, and lead to potentials F{r) and U{r) which are short-ranged [i.e., fall off at a rate 
0(l/r^) or faster as r ^ oo]. In practice, choosing gi{r) = 1 or gi{r) = (r^ + a?)/r'^ appears to lead to well-behaved 
potentials; some experimentation may be needed to make useful choices. 



V. SCALAR RADIATION 

For scalar radiation, s — 0, the functions i?, X, and x have the following relationship: 

X 



go{r)R = X 



A good choice is go{r) = 1. The functions a and (3 [cf. Eq. ( |3.8| )] are then given by 

a = 1 , /? = . 

From this, it follows that 

77=1, 

oi^i = , 
oUi = Vt , 
rA 



qG 



The potentials oF{r) and oU{r) are given by substituting Eqs. ( 5.3a ) - ( 5.3d ) into Eq. (O). For large 



oU{r) = -Lo^ + + ^— '- + 0{l/r^) 



(5.1) 

(5.2) 

(5.3a) 
(5.3b) 
(5.3c) 

(5.3d) 
(5.4) 



clearly, o^'('') = for all r. Hence, the potentials are short-ranged. When a = 0, QU{r) reduces to —lo"^ + V^^{r, s — 0). 
Thus, it reduces to the Regge- Wheeler equation in the Schwarzschild limit, as it was supposed to. 
The asymptotic solutions to the s = GSN equation are simple plane waves: 



X"{r) 



-xpr 



r ^ r-f- 



^°"'Po(r-)e*"'^ -|-^'"Po(r)e- 



The function 



Cn 



Lur (ujr)^ {ujr)^ 



(5.5) 
(5.6) 

(5.7) 



allows us to more accurately describe the behavior of X^'°° near infinity. This is useful both to improve numerical 
computations and to derive certain relations between the amplitudes of the Teukolsky solution and the GSN solution. 
The first three coefficients have the values 



Bo 
Co 



-- (A + 2amuj) , 

-— [A^ — a (2 — Aaniut) — 4 [amoj — iMui — amuj {amoj -t- 2iA/w)]] , 
8 

i 
6 



So (A - 6 + 2amuj + SiMu) ~ 4 {Mujf - 2AoMlu (A - 6) 



(auj)'^ (a — 1 + + 2amcij) 



(5.8) 
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VI. ELECTROMAGNETIC RADIATION 



For electromagnetic radiation, s — —1, the functions R, X, and x exhibit the foUowing relationships: 



X = aR + f3R' 



90 



giir)R 



From this, we can read off 



9i 



rgi igiK 



(6.1) 
(6.2) 
(6.3) 



A 



A useful choice for g^ and gi is 



(r^ + a^) 



9o{r) 



giir) = 1 



The function rj that follows from these choices is 

r] ~ Co + Ci/r + C2/r'^ + Cs/r^ + c^/r^ 



(6.4) 

(6.5) 
(6.6) 



where 



Co 



A 



Using this rj and 



Mr) 



c\ = —2iam , 
C2 = a2(l-2A) , 
C3 = — 2a^(Af + iam) , 
C4 = a*(l - A) . 

rA {r - M) 



(r2 + a2)2 



(6.7) 



(6.8) 



it is straightforward to construct the functions _i-F'i(r), -iU\{r), -iF{r), and -iU{r). The results are rather com- 
plicated and are not given here. When r is large, 



2iam 2a\2imMX- a(2m'^ + 2X'^ - X)] , 



Xr^ ^2^3 
2 A2 + 2amw(A + l) 2 [m2(2m2 - A) + M(A3 + 2amwA)] 



Ar2 



A2r3 



+ 0{\/r^) 



(6.9) 



Both _iF(r) and _i?7(r) are short ranged. When a = 0, _ii^(r) = and _iC/(r) = — tj^ + Vflvy(^, s = ^ !)■ 

Th e solutions X^'°" are, in the limits r — > r_(- and r — > oo, essentially identical to those given in Eqs. (5.5) and 
(5.6); one need only change the subscript on the P function to —1. The corresponding coefficients in P_i(r) are 

A-i ~ — - (A + 2amuj) , 



B-i = -- [A^ - a (2 - Aamuj) - Aauj [m - 2imMui + auj {2 ~ m^)]] 

!_i (A - 6 + 2amuj + 8iMuj) + 2y 
- (aw)2 [A - 5 + SiMuj + 2amLu + m^{X + 2) /A] ] 



C_i = - - [6-1 (A - 6 + 2amuj + 8iMuj) + 2A-i [Mlj (5 - A) + iauj {2au; + m/X)] 



(6.10) 
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VII. CONCLUSION 



In this paper, I have shown how the Sasaki-Nakamura short-ranged equation for gravitational perturbations to a 
rotating black hole may be generalized to arbitrary spin weight radiation. Of course, there is no particular physical 
motivation for choosing radiation spin of magnitude greater than 2; this approach is taken simply so that one can 
write down a single rule which encompasses all physically interesting radiation fields, much as the Teukolsky equation 
itself encompasses all spin weights. Efficient numerical computation of Teukolsky equation solutions now reduces to a 
simple recipe. First, following the analysis in Sec. IV, dev elop the potentials needed in the GSN equation, Eq. ([f.^). 
Examples are given for s = and s = —1. Integrate Eq. ( |4.S| ) for th e G SN solution X. Transform to the variable x 
using Eq. (4.6). Then construct the Teukolsky solution R using Eq. (4.1). 

One application of these results may be to extend the mode sum regularization scheme described in Ref. to 
self forces computed in Kerr spacetimes. Calculations that employ scalar or electromagnetic charges and fields are 
generally simpler than the gravitational self force calculations, which are of great interest for researchers studying 
gravitational-wave sources. The electromagnetic perturbation equation given in Sec. Vl may be of astrophysical 
interest, particularly when coupled to an appropriate source. 
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